Making a field effect transistor on a single graphene nanoribbon by selective doping 
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Using first-principle electronic structure calculations, we show a metal-semiconductor transition 
of a metallic graphene nanoribbon with zigzag edges induced by substitutional doping of Nitrogen 
or Boron atoms at the edges. A field effect transistor consisting of a metal-semiconductor-metal 
junction can then be constructed by selective doping of the ribbon edges. The current-voltage 
characteristics of such a prototype device is determined by the first-principle quantum transport 
calculations. 



Graphene nanoribbons (GNRs) have attracted inten- 
sive interest because of their unique electronic proper- 
ties and vast potential for device applications [l|, 0, H, 
i, i, i, 0, i, M, E [O, [ll- In particular, the GNR- 
based devices could behave like molecule devices, such as 
those based on carbon nanotubes (CNTs)[ll,[l3,[lB|, but 
with some inherent advantages, including more straight- 
forward fabrication processes by using lithography tech- 
nique and better control of crystallographic orientation in 
constructing device junctions,^, St, A JJj . Different from 
CNTs, the existence of edge structures endows GNR with 
some novel physical and chemical properties, such as the 
high edge reactivity and unique edge states around 
the Fermi level[H. These may offer key advantages in re- 
alizing various electronics applications via edge chemical 
functionalization, such as doping. 

It is well known that Nitrogen (N) and Boron (B)atoms 
are typical substitutional dopants in carbon materials 
(like CNTsJJi]), and their binding with the C atom is co- 
valent and quite strong, comparable to that of host C-C 
bond. The incorporation of N or B atoms into the carbon 
materials will influence the electronic and transport prop- 
erties of the C host by introducing extra carries and/or 
new scattering centers [l^. In this letter, we theoretically 
show that a metal-semiconductor transition (MST) can 
be induced in an "armchair" GNR (with zigzag edges) Q 
by substitutional doping of N or B atoms on the edges. 
Based on this finding, we propose a field effect transis- 
tor (FET) made from an single armchair GNR via se- 
lective edge functionalization (doping), and demonstrate 
that the characteristics of such a FET is comparable with 
that of the CNT-FETs. 

Our electronic structure calculations are performed us- 
ing the Vienna Ab-initio Simulation Package [l9l |. which 
implements the formalism of plane wave ultrasoft pscu- 
dopotential based on density functional theory (DFT) 
within local density approximation (LDA). The plane 
wave cut-off energy is set as 350 eV. Structural optimiza- 
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tion was first carried out on all doped systems until the 
residual forces on all ions were converged to below 0.01 
eV/A. The quantum transport calculations were per- 
formed using the Atomistix ToolKit2.0 package [2^, [HI, 
which implements DFT-based real-space, nonequilibrium 
Green's function formalism. The mesh cutoff of carbon 
atom is chosen as 100 Ry to achieve the balance between 
calculation efficiency and accuracy. 

A (2, 2) armchair GNR with H termination (Fig. 
la) is chosen as a model system to study the effect of 
substitutional doping of N and B atoms in armchair 
GNRs. Herein, the "armchair" GNRs are denoted us- 
ing a nomenclature in analogy to armchair carbon nan- 
otubes that would unfold into corresponding ribbons with 
zigzag edges'. Our calculations show that the zero- 
temperature ground state of armchair GNRs is spin- 
polarized in agreement with previous calculations Q. The 
energy of the spin-polarized state is only 20 meV per 
edge atom lower than the spin-unpolarized state. How- 
ever, the spin-polarized state would become unstable 
with respect to the spin-unpolarized state in the presence 
of a ballistic current through the GNRs[l3|- Moreover, 
the magnetization was shown theoretically to be forbid- 
den in pure ID and 2D systems at finite temperatures 
(2^ . while most transistors work at finite temperature 
(room temperature). Therefore, we will only consider 
the spin-unpolarized state of armchair GNRs for our in- 
vestigation of GNR-based devices. In this case, pristine 
armchair (2,2) GNR is metallic with partially fljit bands 
at the Fermi energy (the so-called "edge states" [l] ) local- 
ized at the ribbon edges (Fig. lb). 

Doping is achieved in the supercell made of the 4 unit 
cells by substituting a N or B atom for a C atom in 
the GNRs. Four different doping sites are considered as 
shown in Fig. la, and the corresponding total energies 
are calculated to determine the most energy-favorable 
site. For N doping, the calculated substitution energy of 
site 1 is much lower than those of site 2 (by 1.07 eV), 
site 3 (by 1.00 eV) and site 4 (by 1.32 eV). While for the 
case of B doping, the corresponding energy differences 
are 0.69 eV, 0.60 eV and 0.97 eV, respectively. This 
clearly indicates that the edge (site 1) of GNR is the 
most energetically favorable site for N or B substitution. 
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FIG. 1: (Color online), (a) Atomic structure of the armchair 
(2,2) GNR, where the arrow shows the periodic direction. The 
edges are terminated by the H atoms (denoted by small white 
spheres). Four diflterent substitutional sites are considered, 
(b) Band structure of pure (2,2) GNR. The Fermi energy is set 
to zero and the "edge states" are indicated by the arrows, (c) 
Band structure of the N-doped (2,2) GNR. (d) Band structure 
of the B-doped (2,2) GNR. The partial charge density of two 
bands labelled by the two arrows of the N-doped and B-doped 
GNRs are shown in (e) and (f), respectively. The scale bar is 
in units of e A~"^. 



Furthermore, it is found that the substitution of N or B 
atoms for C atoms does not affect the stability of over- 
all configuration, consistent with previous experimental 
results [231 . The local structural distortion induced by B- 
doping is more pronounced than N-doping, like the case 
in CNTs JJi], which can be related to the atomic radius 
difference between N (or B) and C atom. 

Considering the fact that both edges of a GNR are 
identically active for doping and it is difficult in practice 
to realize selective doping merely on one edge while keep- 
ing another edge unchanged, we will focus our study on 
the cases where two carbon atoms (one per edge) in the 
supercell are simultaneously substituted by two N or B 
atoms. The two substitutional edge sites are determined 
by minimizing the total energy of the system, which are 
indicated by blue atoms in Fig. la. Figs. Ic and Id 
show the band structures of the N-doped and B-doped 
(2,2) GNRs, respectively. From the partial charge den- 
sity analysis (shown in Figs, le and If), we find that 
the energy states (indicated by two arrows in Figs. Ic 
and Id) near the Fermi energy are mainly localized at 
two zigzag edges, and are derived from the original edge 
states. Most importantly, the substitutional doping of 
the N (or B) atoms have removed the degeneracy of the 
eigenvalues at X point and thus open an energy gap. 
The N (B) substitution consequently changes the origi- 
nal band filling [the two states denoted by the arrows are 
now fully occupied (unoccupied)] and eventually induces 
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FIG. 2: (Color online), (a) The dependence of the band gap 
on the linear doping concentration for the N- and B-doped 
(2,2) GNR. (b) The dependence of the band gap on the GNR 
width with the linear doping concentration of 0.1024 A~^. 



a transition of the armchair GNR from metallic to semi- 
conducting. This phenomenon is rather interesting and 
somewhat unexpected since impurity doping, in general, 
results in a transition of semiconducting to metallic. We 
have examined the cases for different substitutional sites 
of two N (B) atoms in GNRs and observed similar MST 
in the system. It should be noted that our above calcula- 
tions correspond to a uniform impurity distribution since 
the conventional periodic boundary condition is adopted 
and there is only one impurity at each edge in the unit 
cell. In order to clarify the effect of random substitution, 
we have also studied electronic structure of GNRs using 
a larger supercell containing two impurities at each edge, 
where the "nonpcriodic" substitution could be partially 
considered with different configurations of the impurities. 
It is found that compared with the periodic substitution 
mentioned above, such "nonperiodic" substitution is en- 
ergetically unfavorable, and importantly, does not change 
the MST of GNRs in essence. The above results show 
that the edge doping-induced transition could be general 
in metallic armchair GNRs. 

Next we will investigate the effect of the doping con- 
centration and nanoribbon width on the electronic prop- 
erties of the doped GNRs. Herein we define the linear 
doping concentration as: ni =Ndopa.nt/ L, where TVdopant 
is the number of dopants per supercell and L is the length 
of the supercell along the nanoribbon. Fig. 2a shows the 
band gaps of N-doped and B-doped (2,2) GNRs as a func- 
tion of the linear doping concentration ni . Note that the 
configuration we studied at each doping concentration is 
always chosen as the energetically most favorable one. 
Interestingly, it is found that the band gap first increases 
with increasing n/ until it reaches the maximum (0.77 eV 
and 0.82 eV for N-doping and B-doping, respectively) at 
ni of 0.1365 A~^, and then decreases with further in- 
creasing ni . The GNRs with B substitution have slightly 
larger band gaps than those with N substitution. In addi- 
tion, the band gap of the GNR decreases with increasing 
ribbon width and eventually diminishes when the width 
is too big(Fig. 2b). 

Among many challenges for the use of graphene or 
nanoribbon in FETs, an important and practical issue 
is to fabricate semiconducting channel with large enough 
band gap, which is crucial for effectively reducing the 
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FIG. 3: (Color online), (a) The schematic structure of the 
field effect transistor (FET) made from a single (2,2) GNR. 
The semiconducting channel is obtained by edge doping of N 
in a finite-length region (the center region), (b) Simulated I — 
1/gato curves of N-doped GNR-FETs under Vtiaa = 0.01 V. The 
channel length is 8.54 nm and the linear doping concentration 
is 0.1365 A~^. (c) The dependence of the subthreshold swing 
S (blue line) and the ON/OFF current ratio (red line) on the 
channel length L. 

leakage current and improving the critical performance 
parameters such as ON/OFF current ratio. However, un- 
til recently it is still very hard to experimentally fabricate 
semiconducting GNRs with the energy gap larger than 
0.2 eV d, The doping-induced MST we report here 
may be used to provide another way to fabricate tran- 
sistor semiconducting channels in future graphene based 
functional devices. Below we will demonstrate the char- 
acteristics and performance parameters of the N-doped 
GNR-FET by first-principle quantum transport calcula- 
tions. It should be noted that both the electrodes and the 
conduction channel are integrated on a single armchair 
GNR in such a device (Fig. 3a). Such a linear config- 
uration can also be advantageous to increase the device 
density in an electronic circuit as well as to simplify the 
fabrication process. 

Fig. 3b shows the typical /-Vgatc curves for the N- 
doped GNR-FET (with the channel length of 8.54 nm) 
at the bias voltage Vbias=0.01 V. In the voltage win- 
dow -0.7 V--0.7 V, the doped FET exhibits ambipolar 
characteristics with the on-current {Ion) of 1 /iA. The 
minimum leakage current is limited to a rather small 
value (-- 1.2x10"^ /iA), and a high ON/OFF current 
ratio {lon/IoS > 2000) is achieved in such N-doped 
GNR-FETs. The large ON/OFF current ratio mani- 



fests the "perfect" atomic interface between the metal- 
semiconductor GNR junctions with a minimum con- 
tact resistance. It increases the possibility of experi- 
mental operation between ON and OFF states. More- 
over, such a device exhibits an excellent "theoretical" 
switching characteristics with a subthreshold swing S — 
ln(10)[dygato/d(ln/)] - 40 mV/decade. This is com- 
parable to that of high performance CNT-FET (60-80 
m V/decade) [H [isj and reaches the theoretical limit of 
S (~ 60 mV/decade) for Si-based FET at room temper- 
ature [2^ . The switching mechanism here can be under- 
stood in terms of a semiclassical band bending mecha- 
nism @. 

We further studied the relationship between the de- 
vice performance and the channel length by calculating 
I-Vgato curve of N-doped GNR-FETs as a function of 
the doped channel length from 0.49 nm to 8.54 nm while 
keeping the bias voltage Vbias at 0.01 V. As shown in Fig. 
3c, the subthreshold swing S of these doped GNR-FETs 
decreases and the ON/OFF current ratio increases expo- 
nentially. Our calculations show that in order to obtain 
good device performance with small S value (e. g., below 
100 mV/decade) and high ON/OFF current ratio (e. g., 
above 100), the doped channel length needs to be longer 
than 5 nm. The minimum leakage current of those FETs 
with the doped channels shorter than this critical length 
will be greatly enhanced by direct tunneling, which low- 
ers the device performance. 

In conclusion, using first-principle calculations, we 
have studied the electronic and transport properties of 
armchair GNRs with substitutional edge doping of N or 
B atoms. It is found that the edge doping will greatly 
modify the band structure (especially the "edge states" ) 
of the system, and induce a metal-semiconductor tran- 
sition. The band gap of the doped GNR exhibits a 
strong dependence on both the linear doping concentra- 
tion and the nanoribbon width. It is demonstrated that 
electronic devices, such as FETs, could be integrated on 
a single GNR by selective edge doping/chemical func- 
tionalization, which avoids the need to connect/integrate 
the GNRs with different orientations. Simulated /-V^ate 
curves indicate that such FETs exhibit ambipolar ON- 
OFF characteristics with excellent performance parame- 
ters. 
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